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We study the effects of higher order transversal modes in a model of a singly- 
resonant OPO, using both numerical solutions and mode expansions including up 
to two radial modes. The numerical and two-mode solutions predict lower 
threshold and higher conversion than the single-mode solution at negative 
dispersion. Relative power in the zero order radial mode ranges from about 88% at 
positive and small negative dispersion to 48% at larger negative dispersion, with 
most of the higher mode content in the first mode, and less than 2% in higher 

modes. 
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Continuous-wave operation of singly-resonant optical parametric oscillators (SRO) 
has recently been demonstrated for various non-linear materials, including KTP [1- 
3] and periodically-poled LiNb03 (PPLN) [4-6]. Significant pump depletion was 

achieved, up to 93% in [6]. In these experiments, with zero beam walk-off and 
confocal parameters of pump and signal comparable to the crystal length, 
diffractive effects are known to be important [7,8]. Theoretical studies including 
diffractive effects for cw SROs [7,8] have been published for the case of small 
pump depletion, and give the dependence of the threshold on the focusing 
geometry. A treatment of a low-loss SRO for arbitrary pump depletion in the 
plane-wave approximation was given in [9]. A numerical model including 
diffraction and pump depletion for a nano-second pulsed SRO with significant 
birefrigent walk-off was presented in [10]. To the best of our knowledge, no study 
of diffractive effects for arbitrary pump depletion in a cw SRO without beam walk- 
off has been published so far. In particular, it is not known how the transverse 
mode content of the idler beam and the amount of pump depletion depend upon 
pump power and dispersion. In this letter, we present such results for an SRO with 
low cavity losses [2]. We find that the idler beam contains higher transverse modes 
with appreciable amplitudes, and that these amplitudes depend only weakly on the 
amount of pump depletion, but strongly on the amount of dispersion. 

We use the coupled wave equations for signal, idler and pump fields Ap 

(F=S,I,P) in the form: 



A S , Z = G/2ks)(A S ,xx + As, yy ) + j K S exp(jAz)A P Ai 



* 



(La) 



Al, z = G/2kl)(Ai 5XX + Ai 5yy ) + j ki expGAz)A P A S 



(Lb) 



Ap, z = G/2kp)(Ap 5XX + Ap 5yy ) + j k P exp(-jAz)A S A! . 



(l.c) 



where Kp=cop d/(np c) and d is the effective nonlinear constant. The subscripts 
x,y and z denote partial derivatives with respect to these spatial coordinates and, 



for a periodically-poled crystal, A=kp-k§-kj - 271/1 is the residual dispersion wave 

vector ( 1 is the period of the poling). In the limit of small cavity loss «s «!> tne 

signal field has almost constant power. For a Gaussian pump input, the lowest 
transverse signal mode has largest gain, and therefore reaches threshold first. We 
found that the threshold for the first higher radial signal mode is about four times 
higher than for the lowest mode. Our model therefore approximates the signal field 
by a Gaussian with constant power . In this case, the relationship between pump 
input power and idler output power is linear, to zeroth order in the cavity loss «s • 

From photon conservation, the first order increment in signal power is then found 
to be proportional to the idler output and hence to the pump input. Equating this to 
the round-trip cavity loss determines the pump input power as a function of signal 
power [2] . In the limit of small signal power, or small pump depletion, this gives 
the threshold pump power Pthr- 

In order to understand the influence of dispersion and pump power on the mode 
content of the idler and pump, it is instructive to consider a mode expansion of the 
fields Ap into radially symmetric modes Bp n [11,12]. For a zero order signal 

mode with constant coefficient, the following coupled mode equations result [12]: 

Ap(x,y,z) = £apn(z) Bp n (x,y,z) , 
n 

dai n /dz = Z c nm ap m aso* , (2.a) 
m 

dap m /dz = Z c nm * ai n aso • (2.b) 
n 



If pump, signal and idler modes have equal confocal parameter and waist 



location, the coupling coefficient c nm takes the form 



Cnm = dmn exp(j<J>)/V l+(z/z ) 2 , (3) 

= zA + (2n+l-2m)arctan(z/zo) = ( A'+l+2n-2m) (z/zo), 
dio = (1-Xp/Xi) doo , di i = (kp/Xi) doo , doi = . 

Here zq is the Raleigh range, A' = zq A, and the coefficient doo is independent of 
z. The approximate form of the phase angle given above is valid for z/zo « 1- 

At threshold and with a zero order pump at input (m=0), the effect of the z- 
dependent phase factor is to favor the n=0 idler mode for A' = -1 and the n=l mode 
for A' = -3. Lowest threshold will occur for -3 < A' < -1. The amplitude of the zero 
idler mode (n=0) will be largest, the amplitude of the first mode (n=l) will be 
smaller, while the amplitudes of higher modes (n>l) will be suppressed by the 
larger phase mismatch. 

Above threshold, the mode content of the idler will be modified by the higher 
pump modes (m>0) generated by pump depletion. However, since d()m=0, such 

pump modes are generated only by higher signal modes (n>0). Unless A' assumes 
large negative values, these have lower amplitude. Furthermore, the coefficients 
d nm , n>0, are small. Thus, we expect the idler mode content to depend only 

weakly on the pump power. 

Our model calculations involve a numerical integration of the coupled wave 
equations 1 .b and 1 .c, and of the coupled mode equations 2. The integration of 
equations 2 uses a second order Runge Kutta algorithm, for the case of zero order 
modes only (m,n=0; single mode approximation SMA) and for both the zero and 
first order radial modes (m,n=0,l; two mode approximation TMA). 

The integration of 1 .b and 1 .c uses a split-step algorithm, in which an explicit 
half-step of parametric conversion is followed by a full step of diffraction and 
implicit half-step of conversion. The diffraction step uses the Cayley-transform 
[13] of a discretized version of the radial Laplacian, which conserves photon 



number exactly. The program was validated by propagating Gaussian modes 
without conversion, and beams with negligible diffraction and full conversion. 
Other numerical results were validated by reducing the radial and longitudinal 
step-sizes by a factor 0.5, giving agreement of better than 0.1% in beam 
amplitudes. 

The models refer to a traveling wave cavity, zero walk-off, equal confocal 
parameter for pump and signal, crystal length L equal to the confocal parameter b 
= 2zq and waist locations at the center of the crystal. Wavelengths are 1 .064 [im 

for the pump and 3.25 [im for the idler. Pump power P and dispersion A are 
reported as dimensionless parameters P' and A'=zo A: 

P = cx S P0 P'. 

P0 = eo c np ns nj c-^ /( 32 p^ cos & z kp ). 

With this unit of power, the threshold parameter P^r' = Pthr/( a SPo) is related to 
the threshold function h of [8] by: 

P thr ' =0.25 (z /L)(l+Xs/Xp)/h . 

For the SMA, the threshold parameter P^r' is given in terms of the Boyd- 
Kleinman gain reduction factor h m [7] as: 

Pthr' = 0.5 (zo*L)(A,i ^P 2 )/hm . 

Figure 1 shows the threshold parameter P^r' as a function of the dispersion 

parameter A'. The SMA gives minimal threshold at A'=-0.86, in good agreement 
with the discussion above. The TMA and the numerical model give minimal 



threshold parameter P' for A'=-1.30. The TMA agrees very well with the numerical 
model for the range of A' considered. The SMA agrees with these for A'>-0.4, but 
gives higher threshold parameter for smaller values of A'. The minimal threshold 
parameter is P' = 1.19 for the numerical model and the TMA, and 1.46 for the 
SMA. The corresponding values for the gain reduction factors h and h m are h = 

0.26 and h m = 0.778. This agrees with the published values in [8] and [7]. 

To determine the actual threshold, the wavelength dependence of both the 
dispersion parameter A' and the resonator loss ocs needs to be considered. The 

signal wavelength will adjust to a value which gives lowest threshold Pthr- 

ocsPthr'PO • If a S i s > within a certain range, independent of wavelength, minimal 

threshold occurs at a wavelength giving the optimal value of A' = -1.30. If «s has a 

minimum at some other wave length, minimal threshold will occur at 
correspondingly different values of A'. 

Wavelength acceptance bandwidths for variations in A' have been given in [4] 
for a PPNL crystal of length L=lcm. For the d^3 coefficient and a shift in A' of 1.0 

above and below the optimal value, the acceptance bandwidth is 7 nm. According 
to the plane-wave approximation (PWA)[9], such a shift leads to an increase in 
threshold power by a factor of 1.41. Our numerical model predicts a smaller 
increase, by a factor of about 1.20. 

Figure 2 shows the maximal pump depletion as a function of dispersion. The 
SMA predicts 100% depletion at A'=-0.86 and a pump power of 2.47 times above 
threshold. This is the same pump power ratio as that predicted by the PWA [9] 
Both the numerical model and the TMA predict maximal depletion of about 95% at 
A'=-1.30 and a pump power of 2.67 above threshold. This is in good agreement 
with the experiment reported in [6], which found maximal depletion of about 93% 
between 2.6 and 3.7 times above threshold. The optimal values of the dispersion 
parameter A' are the same as those found for minimal threshold. As above, 
numerical model and TMA agree very well, while the SMA predicts lower 
conversion and higher pump requirements for A'<-0.8. For a shift in A' by 1.0 



below and above optimum, our numerical model predicts a reduction of maximal 
pump depletion to 81% and 75%, while the SMA gives 69% and 68%. The PWA 
gives 66% in both cases. 

These reductions in threshold and increases in pump depletion, as compared to 
the PWA and SMA, are a consequence of the variable mode structure of the idler 
beam, as discussed above. As A' is shifted away from the optimal value, any 
reduction in gain for one of the modes is partially offset by an increase in gain for 
other modes, leading to a strong dependence of the idler mode content on A'. 

Figure 3 shows the relative power in the idler modes n=0,l and n>l at threshold 
as functions of A'. As predicted, the n=0 mode dominates for small negative and 
positive values of A', while the n=l mode becomes increasingly stronger for larger 
negative values, and dominates for A'<-2.4. For the optimal value A'=-1.3, the 
mode n=0 contains 77% of the idler power, the mode n=l contains 22%, and 
modes n>l contain 1%. If the focal point z=z'zq and the Raleigh range zj=zj'zq 

of the idler modes are adjusted, to minimize the relative power in modes n>l, such 
modes contain no more than 2% of the total power. For example, at A'=-2.3, for 
modes with focus at z'=1.0 and zj'=1.0, the relative power in modes n=0,l and n>l 

is 51%, 44% and 5%, while for the adjusted modes with focus at z'=1.04 and 
zj'=0.74, the relative power is 65%, 34% and 1%. At A'=-1.3, the adjusted modes 

have focus at z'=1.04, Raleigh range zj'=0.94 and contain 80% for n=0, 20% for 

n=l and less than 1% for n>l. These values are very close to the values obtained 
for the original set of modes.. 

Figure 4 finally shows the idler mode content above threshold as a function of 
pump power, for adjusted modes, at A'=-2.3, -1.3 and -0.3. Fig. 4.a shows the 
relative power in the mode n=l, fig. 4.b shows the relative phase of modes n=0 and 
n=l in degrees, and fig. 4.c shows the location of the mode waist. The relative 
power in the mode n=l is seen to be almost independent of pump power, but 
depends strongly on dispersion. This strong dependence of higher mode amplitudes 
on the dispersive phase mismatch was also observed in the experimental and 
numerical studies reported in [10]. The relative phase, which determines the beam 



radius in the far-field region, vanishes at threshold and increases in absolute value 
with pump power. This leads to an increase of the far-field beam radius with pump 
power. The location of the waist shifts towards the entry face of the crystal (z'=0), 
as the pump power increases. The Raleigh range of the adjusted modes is almost 
independent of pump power, but strongly dependent upon dispersion. It ranges 
from zj'= 0.74 at threshold to 0.72 at maximal depletion for A'=-2.3, from 0.94 to 

0.93 for A'=-1.3 and from 1.24 to 1.26 at A'=-0.3. 

In summary, the models studied here demonstrate that the idler output of a cw, 
low loss SRO with collinear beams contains higher transverse modes with 
appreciable amplitude, with most of the higher mode content in the first radial 
mode. The mode content depends strongly on dispersion, but only weakly on pump 
power. Furthermore, the variable mode content of the idler leads to lower threshold 
and in higher conversion efficiency for larger negative values of the dispersion 
parameter A' than those predicted by the single-mode, or plane-wave 
approximations. The price paid for this is a higher content of radial modes n>0 in 
the idler output. 

It is a pleasure to acknowledge Walt Bosenberg at Lightwave Electronics for 
discussions about the subject of this paper, and for making preprints available. 
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Figure Captions 

Figure 1 : Threshold pump power as a function of the dispersion parameter A'=zqA, 
for the numerical, two-mode and single-mode models. 

Figure 2: Pump depletion as a function of the dispersion parameter A' at maximal 
conversion, for the numerical, two-mode and single-mode models. 

Figure 3: Relative power in the idler modes n=0, n=l and n>l at threshold, as a 
function of the dispersion parameter A' 

Figure 4: Mode content of the idler beam as function of normalized pump power 
P'=P/Pthr> f° r A'=-2.3, -1.3 and -0.3. Mode parameters are adjusted for minimal 

power in modes n>l. Fig. 4.a shows the relative power in mode n=l, fig. 4.b shows 
the relative phase of mode n=l, and fig. 4.c shows the location z'=z/zq of the 

adjusted mode waist. 
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